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ABSTRACT 

We  present  three  complete  seasons  and  two  half-seasons  of  Sloan  Digital  Sky  Survey  (SDSS)  r-band  photometry 
of  the  gravitationally  lensed  quasar  SBS  0909+532  from  the  U.S.  Naval  Observatory,  as  well  as  two  seasons  each 
of  SDSS  g-band  and  r-band  monitoring  from  the  Liverpool  Robotic  Telescope.  Using  Monte  Carlo  simulations  to 
simultaneously  measure  the  system’s  time  delay  and  model  the  r-band  microlensing  variability,  we  confirm  and 
significantly  refine  the  precision  of  the  system’s  time  delay  to  A t^B  —  50+24  days,  where  the  stated  uncertainties 
represent  the  bounds  of  the  formal  1  a  confidence  interval.  There  may  be  a  conflict  between  the  time  delay 
measurement  and  a  lens  consisting  of  a  single  galaxy.  While  models  based  on  the  Hubble  Space  Telescope 
astrometry  and  a  relatively  compact  stellar  distribution  can  reproduce  the  observed  delay,  the  models  have  somewhat 
less  dark  matter  than  we  would  typically  expect.  We  also  carry  out  a  joint  analysis  of  the  microlensing  variability 
in  the  r  and  g  bands  to  constrain  the  size  of  the  quasar’s  continuum  source  at  these  wavelengths,  obtaining 
logjOyr/cmjtcosi'/O.S]1^2}  =  15.3  ±  0.3  and  log{(rJig/cm)[cos  f  /0.5] 1/2}  =  14.8  ±  0.9,  respectively.  Our  current 
results  do  not  formally  constrain  the  temperature  profile  of  the  accretion  disk  but  are  consistent  with  the  expectations 
of  standard  thin  disk  theory. 

Key  words:  accretion,  accretion  disks  -  gravitational  lensing:  micro  -  gravitational  lensing:  strong  -  quasars: 
individual  (SBS  0909+532) 


1.  INTRODUCTION 

Much  of  the  standard  picture  for  the  detailed  structure  of 
accretion  disks  surrounding  supermassive  black  holes  in  active 
galactic  nuclei  is  based  on  theoretical  models  rather  than  ob¬ 
servational  measurements,  because  these  very  compact  regions 
cannot  be  resolved  with  existing  telescopes.  Yet,  for  gravita¬ 
tionally  lensed  quasars,  the  relative  motions  of  the  observer, 
the  background  source,  the  foreground  lens  galaxy,  and  its  stars 
cause  uncorrelated  variations  in  the  source  magnification  as  a 
function  of  time  and  wavelength  which  depend  on  the  projected 
area  of  the  continuum  source.  By  analyzing  these  microlens¬ 
ing  brightness  fluctuations  with  numerical  simulations,  one  can 
measure  the  continuum  source  size,  permitting  invaluable  ob¬ 
servational  tests  of  theoretical  models  of  accretion  disk  structure 
(e.g.,  Eigenbrod  et  al.  2008;  Anguita  et  al.  2008;  Poindexter  et  al. 
2008)  and  orientation  (Poindexter  &  Kochanek  2010).  Such  tests 
in  the  literature  typically  fall  into  two  categories.  One  type  of 
study  employs  single-epoch  multi-band  photometry  of  lensed 
quasars  in  which  flux  ratios  of  the  images  exhibit  deviations 
from  the  predictions  of  macroscopic  lens  models  or  significant 
wavelength  dependence  (e.g.,  Pooley  et  al.  2006;  Bate  et  al. 
2008;  Blackburne  et  al.  2011;  Mosquera  et  al.  2011;  Motta  et  al. 
2012).  The  alternative  method,  described  in  detail  in  Kochanek 
(2004),  analyzes  the  time  variability  of  the  quasar’s  flux  ratio  and 
requires  monitoring  a  quasar  over  a  significant  period  of  time. 
Although  such  investigations  can  be  observationally  and  com¬ 
putationally  challenging,  they  have  the  advantage  that  it  is  not 


necessary  to  assume  a  value  for  the  mass  of  the  microlenses  or 
an  extinction  law  for  the  lens  galaxy,  nor  are  the  results  highly 
sensitive  to  assumed  priors.  Studies  employing  both  methods 
of  analysis  have  revealed  accretion  disk  temperature  profiles 
in  agreement  with  the  simple  thin-disk  model  of  Shakura  & 
Sunyaev  (1973),  although  the  individual  observed  disk  sizes 
tend  to  be  larger  than  those  predicted  by  the  theoretical  model 
(e.g.,  Pooley  et  al.  2007;  Morgan  et  al.  2010;  Blackburne  et  al. 
2011). 

SBS  0909+532  (hereafter  SBS  0909;  aj2000  =  09h13m0L05, 

5j7000  =  +52c*59m28)83)  is  a  doubly-imaged  quasar  lens  sys¬ 
tem  in  which  the  background  quasar  has  redshift  zs  =  1 .377  and 
the  foreground  early-type  lens  galaxy  has  redshift  zi  =  0.830 
(Kochanek  et  al.  1997;  Oscoz  et  al.  1997;  Lubin  et  al.  2000). 
SBS  0909  is  a  somewhat  challenging  system  to  study  be¬ 
cause  discrepant  results  in  the  literature  expose  significant 
uncertainties  about  some  of  its  fundamental  properties  and 
the  nature  of  the  quasar’s  variability.  Most  notably,  Lehar 
et  al.  (2000)  encountered  difficulty  when  attempting  to  use 
the  imfitfits  routine  to  measure  the  lens  galaxy  photometry  in 
//-band  NICMOS  images  from  the  Hubble  Space  Telescope 
(HST),  eventually  settling  on  a  low  surface  brightness  de  Vau- 
couleurs  model  with  effective  radius  repp  =  1"58  ±  O'.' 9  and 
magnitude  H  =  16.75.  In  stark  contrast,  Sluse  et  al.  (2012) 
used  an  iterative  deconvolution  technique  (e.g.,  Magain  et  al. 
1998;  Chantry  &  Magain  2007)  on  the  same  data  to  find  a 
significantly  smaller  (repp  =  (X'54  ±  O'.' 02)  and  less  luminous 
(H  =  19.44  ±0.01)  lens  galaxy,  formally  inconsistent  with 
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the  Lehar  et  al.  (2000)  result.  Additionally,  past  optical  mon¬ 
itoring  observations  have  not  shown  evidence  of  significant 
microlensing  variability  (Ullan  et  al.  2006;  Goicoechea  et  al. 
2008),  although  Mediavilla  et  al.  (201 1)  find  evidence  of  chro¬ 
matic  microlensing  (wavelength-dependent  microlensing  mag¬ 
nification)  through  an  analysis  of  the  differences  between  the 
continuum  and  emission  line  flux  ratios  from  the  quasar’s  ultra¬ 
violet  (UV)-near-infrared  (near-IR)  spectra,  which  can  separate 
microlensing  effects  from  differential  extinction  present  in  the 
lens  galaxy. 

We  have  compiled  a  new  data  set  consisting  of  monitoring 
observations  of  SBS  0909  in  two  optical  bands  from  two 
different  observatories,  the  analysis  of  which  provides  some 
resolution  to  these  discrepant  results  from  the  literature.  Here  we 
will  show  that  in  the  four  years  spanning  2008-2012,  SBS  0909 
has  exhibited  significant  uncorrelated  time  variability  in  the  rest- 
frame  near-UV.  We  will  analyze  the  uncorrelated  variability, 
which  we  attribute  to  microlensing,  to  determine  a  size  for  the 
accretion  disk  in  the  two  different  bands,  and  then  compare  our 
results  to  that  derived  from  the  observed  chromatic  microlensing 
by  Mediavilla  et  al.  (2011).  Since  our  multi-band  data  set  allows 
us  to  constrain  the  size  of  the  accretion  disk  at  two  different 
wavelengths,  we  also  gain  a  glimpse  at  the  temperature  profile 
of  the  accretion  disk,  with  less  reliance  on  priors  than  single¬ 
epoch  microlensing  analyses. 

In  order  to  analyze  the  uncorrelated  variability  of  a  lensed 
quasar  using  the  methods  of  Kochanek  (2004),  we  must  first 
have  accurate  knowledge  of  the  time  delay  between  the  lensed 
images.  In  general,  imprecise  time  delays  can  result  in  signifi¬ 
cant  uncertainties  in  microlensing  analyses,  since  residual  vari¬ 
ability  from  an  improperly  corrected  time  delay  must  be  mod¬ 
eled  as  uncorrelated  microlensing  variability  in  simulated  light 
curves.  In  the  case  of  SBS  0909,  the  lens  mass  model  provides 
very  little  help  with  reducing  the  time  delay  uncertainty  given 
the  very  discrepant  lens  galaxy  photometric  fits  of  Lehar  et  al. 
(2000)  and  Sluse  et  al.  (2012).  Unfortunately,  the  time  delay  of 
SBS  0909  has  been  rather  difficult  to  determine  because  there 
have  been  relatively  few  large-amplitude  (>0. 1  mag),  short- 
duration  extrema  in  the  light  curves  in  the  years  since  its  dis¬ 
covery.  The  most  recent  time  delay  measurement  for  this  quasar, 
A tAB  —  49  ±  6  days  (Goicoechea  et  al.  2008),  still  had  relatively 
significant  uncertainty.  In  our  new  /  -band  data  set  for  SBS  0909, 
the  quasar  images  exhibit  substantial  intrinsic  flux  variation,  but 
this  flux  variation  is  modulated  by  the  uncorrelated  microlensing 
variability  so  that  a  simple  polynomial-based  cross-correlation 
analysis  (e.g.,  Kochanek  et  al.  2006;  Poindexter  et  al.  2007)  fails, 
a  problem  explored  by  Eigenbrod  et  al.  (2005).  Morgan  et  al. 
(2008)  introduced  a  method  of  making  this  problem  tractable  by 
analyzing  the  microlensing  variability  with  the  techniques  from 
Kochanek  (2004)  while  simultaneously  solving  for  the  time  de¬ 
lay.  We  will  apply  this  methodology  to  our  /  -band  data  set  to 
make  a  new,  independent  measurement  of  the  time  delay  for 
SBS  0909  which  provides  strong  evidence  in  favor  of  a  more 
compact  lens  galaxy  photometric  model  (e.g.,  Sluse  et  al.  2012) 
but  at  the  astrometric  position  favored  by  Lehar  et  al.  (2000). 

In  Section  2  we  describe  our  g-  and  /  -band  monitoring  obser¬ 
vations  from  the  United  States  Naval  Observatory — Flagstaff 
and  the  Liverpool  Telescope  and  how  we  combined  the  two  data 
sets.  In  Section  3  we  outline  our  Bayesian  Monte  Carlo  method 
for  simultaneously  estimating  time  delays  and  quasar  structure. 
In  Section  4  we  present  the  results  of  our  analysis  and  compare 
our  findings  to  those  of  Mediavilla  et  al.  (2011).  Throughout 
our  discussion,  we  assume  a  flat  cosmology  with  =  0.3, 


£2a  =  0.7,  and  Hq  =  70kms  1  Mpc  1  (Hinshaw  et  al. 
2009). 

2.  OBSERVATIONAL  DATA 
2.1.  USNO  Monitoring 

We  observed  SBS  0909  regularly  as  part  of  the  United 
States  Naval  Academy /United  States  Naval  Observatory 
(USNA/USNO)  Lensed  Quasar  Monitoring  Program.  Using 
the  1.55  m  Kaj  Strand  Astrometric  Reflector  at  the  USNO — 
Flagstaff  Station,  we  take  three  five-minute  exposures  of  the 
quasar  in  the  Sloan  Digital  Sky  Survey  (SDSS)  /--band  per 
epoch  at  a  cadence  of  two  to  three  nights  per  month,  using 
either  the  2048  x  2048  Tek2K  CCD  camera  (O'.' 3 3  pixel-1)  or 
the  2048  x  4096  EEV  CCD  camera  (O'.' 18  pixel-1).  The  de¬ 
tails  of  our  photometric  analysis  techniques  are  discussed  in 
detail  in  Kochanek  et  al.  (2006).  In  summary,  we  measure  the 
quasar  image  fluxes  relative  to  three  reference  stars,  located  at 
(—12/8,  71"9),  (67" 2,  11/4),  and  (-13/8,  7/1)  with  respect  to 
image  A  of  SBS  0909,  using  a  three-component  elliptical  Gaus¬ 
sian  point-spread  function  (PSF)  model.  We  hold  the  relative 
positions  of  the  quasar  images  fixed  to  those  derived  from  the 
HST  H- band  images  of  SBS  0909  for  the  PSF  fitting  process. 
The  photometric  model  of  the  very  red  zi  =  0.83  lens  galaxy  is 
a  Gaussian  approximation  to  a  de  Vaucouleurs  profile  of  fixed 
effective  radius  and  flux.  We  use  the  effective  radius  derived 
from  the  HST  images  by  Lehar  et  al.  (2000)  and  for  the  flux  we 
use  the  value  which  minimizes  the  total  x  2  in  the  residuals  over 
all  epochs.  We  also  attempted  to  measure  the  quasar  photometry 
using  the  more  compact  and  dimmer  lens  galaxy  photometric 
fit  from  Sluse  et  al.  (2012),  but  the  changes  in  the  quasar  im¬ 
age  fluxes  were  negligible.  Unlike  the  case  of  Q  0957+561  in 
Hainline  et  al.  (2012),  no  color  offset  is  required  between  the  two 
different  detectors  used  for  our  observing  program.  In  Table  1, 
we  present  the  /-band  measurements  of  SBS  0909  components 
A  and  B  from  61  nights  between  2008  March  and  2012  Febru¬ 
ary.  The  images  from  which  our  measurements  are  derived  are 
characterized  by  a  median  stellar  FWHM  (seeing)  of  1/3.  Be¬ 
cause  the  quasar  images  are  closely  spaced  (1/17),  they  are 
blended  in  our  USNO  images,  causing  our  photometric  analysis 
to  break  down  for  seeing  conditions  somewhat  larger  than  the 
image  separation,  so  we  keep  only  epochs  for  which  the  see¬ 
ing  is  better  than  1/6  in  our  analysis.  This  removes  10  epochs 
from  our  USNO  data  set,  which  are  identified  in  Table  1  by  the 
parentheses  surrounding  the  measurements.  We  discarded  an 
additional  six  epochs  not  listed  in  Table  1  due  to  partial  cloud 
cover  and  bright  sky  conditions. 

2.2.  Liverpool  Telescope  Monitoring 

We  also  monitored  SBS  0909  in  the  r  band  with  the  2.0  m 
Liverpool  Robotic  Telescope  (LRT)  independently  of  the  ob¬ 
servations  at  USNO.  Our  /  -band  monitoring  program  used  the 
RATCam  CCD  camera,  providing  a  4 .'6  field  of  view  with  pixel 
scale  0/28  pixel-1,  and  was  carried  out  over  two  different  pe¬ 
riods:  from  2005  January  to  2007  January  (I),  and  from  2010 
October  to  2012  March  (II).  The  measurements  spanning  the 
time  period  2005  January-2006  June  (78  epochs)  have  already 
been  published  in  Goicoechea  et  al.  (2008).  Here  we  add  30 
additional  epochs  of  magnitudes  corresponding  to  the  last  time 
segment  from  monitoring  period  I  (2006  October  to  2007  Jan¬ 
uary)  and  the  first  90  epochs  from  monitoring  period  II  (2010 
October  to  20 1 1  June),  for  a  total  of  1 98  epochs  on  the  Liverpool 
Telescope. 
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Table  1 

SBS  0909+532  Light  Curves 


HJD  -  2450000 
(days) 

Seeing 

(arcsec) 

QSO  A 
(mag) 

QSO  B 
(mag) 

(Stars) 

(mag) 

4554.635 

1.2 

2.475 

± 

0.005 

3.281 

± 

0.006 

0.033 

± 

0.003 

4555.635 

1.2 

2.483 

± 

0.005 

3.279 

± 

0.006 

0.034 

± 

0.003 

4561.686 

1.7 

2.483 

± 

0.005 

3.277 

± 

0.006 

0.031 

± 

0.003 

4570.625 

1.5 

2.490 

± 

0.006 

3.249 

± 

0.009 

-0.087 

± 

0.003 

(4584.638) 

(1.9) 

(2.494 

± 

0.005) 

(3.245 

± 

0.006) 

(0.008 

± 

0.003) 

4596.651 

1.3 

2.493 

± 

0.005 

3.273 

± 

0.007 

-0.016 

± 

0.003 

4613.664 

1.6 

2.529 

± 

0.009 

3.227 

± 

0.015 

-0.061 

± 

0.004 

4617.673 

1.5 

2.515 

± 

0.006 

3.229 

± 

0.007 

0.014 

± 

0.003 

4627.656 

1.3 

2.484 

± 

0.006 

3.278 

± 

0.010 

-0.120 

± 

0.003 

4792.967 

1.0 

2.504 

± 

0.007 

3.188 

± 

0.009 

-0.061 

± 

0.004 

4802.946 

1.2 

2.511 

± 

0.005 

3.168 

± 

0.006 

0.013 

± 

0.003 

4807.895 

1.3 

2.517 

± 

0.006 

3.120 

± 

0.008 

-0.037 

± 

0.003 

4829.824 

1.3 

2.522 

± 

0.005 

3.125 

± 

0.006 

0.016 

± 

0.003 

4833.850 

1.2 

2.469 

± 

0.005 

3.143 

± 

0.006 

0.022 

± 

0.003 

4839.929 

1.3 

2.478 

± 

0.006 

3.131 

± 

0.007 

-0.035 

± 

0.003 

4862.837 

1.4 

2.456 

± 

0.005 

3.171 

± 

0.006 

0.046 

± 

0.003 

4883.892 

1.6 

2.489 

± 

0.008 

3.060 

± 

0.011 

-0.052 

± 

0.003 

4891.796 

1.3 

2.481 

± 

0.006 

3.073 

± 

0.007 

-0.009 

± 

0.003 

4911.756 

1.0 

2.471 

± 

0.005 

3.070 

± 

0.005 

0.039 

± 

0.003 

4942.748 

1.3 

2.453 

± 

0.005 

3.060 

± 

0.006 

0.023 

± 

0.003 

4949.723 

1.5 

2.453 

± 

0.006 

3.046 

± 

0.007 

-0.005 

± 

0.003 

4964.683 

1.2 

2.437 

± 

0.004 

3.043 

± 

0.005 

0.042 

± 

0.002 

5157.914 

1.4 

2.425 

± 

0.005 

2.986 

± 

0.006 

0.021 

± 

0.003 

5163.034 

1.3 

2.430 

± 

0.005 

2.962 

± 

0.005 

0.034 

± 

0.003 

5185.971 

1.4 

2.413 

± 

0.006 

2.927 

± 

0.007 

-0.005 

± 

0.003 

5201.914 

1.3 

2.406 

± 

0.005 

2.920 

± 

0.006 

-0.028 

± 

0.003 

5208.960 

1.2 

2.409 

± 

0.005 

2.894 

± 

0.005 

0.024 

± 

0.003 

(5241.774) 

(1.7) 

(2.385 

± 

0.005) 

(2.861 

± 

0.006) 

(0.006 

± 

0.002) 

5296.650 

1.0 

2.356 

± 

0.004 

2.889 

± 

0.004 

0.059 

± 

0.002 

5311.656 

1.4 

2.352 

± 

0.006 

2.888 

± 

0.007 

-0.074 

± 

0.003 

5324.698 

0.8 

2.357 

± 

0.005 

2.866 

± 

0.006 

0.037 

± 

0.003 

5332.680 

1.1 

2.363 

± 

0.005 

2.896 

± 

0.005 

0.024 

± 

0.003 

(5339.656) 

(1.7) 

(2.374 

± 

0.007) 

(2.875 

± 

0.009) 

(-0.086 

± 

0.003) 

5348.666 

1.0 

2.350 

± 

0.004 

2.905 

± 

0.004 

0.043 

± 

0.002 

5354.659 

1.1 

2.351 

± 

0.004 

2.891 

± 

0.005 

-0.002 

± 

0.002 

5471.991 

1.0 

2.344 

± 

0.004 

2.900 

± 

0.004 

0.028 

± 

0.002 

5481.013 

1.2 

2.350 

± 

0.005 

2.906 

± 

0.006 

0.012 

± 

0.003 

5484.007 

1.2 

2.330 

± 

0.005 

2.940 

± 

0.005 

0.030 

± 

0.003 

5506.953 

1.4 

2.351 

± 

0.005 

2.917 

± 

0.006 

0.008 

± 

0.003 

5519.016 

1.1 

2.354 

± 

0.005 

2.895 

± 

0.005 

0.032 

± 

0.003 

5534.989 

1.2 

2.350 

± 

0.006 

2.866 

± 

0.007 

-0.008 

± 

0.003 

5543.959 

1.4 

2.354 

± 

0.005 

2.878 

± 

0.005 

0.025 

± 

0.003 

5563.891 

1.3 

2.325 

± 

0.004 

2.876 

± 

0.004 

0.025 

± 

0.002 

5590.748 

1.1 

2.315 

± 

0.004 

2.898 

± 

0.004 

0.048 

± 

0.002 

(5604.800) 

(1.9) 

(2.331 

± 

0.005) 

(2.897 

± 

0.006) 

(0.009 

± 

0.003) 

5621.653 

1.4 

2.318 

± 

0.005 

2.891 

± 

0.006 

-0.067 

± 

0.003 

5626.747 

1.6 

2.326 

± 

0.005 

2.901 

± 

0.006 

0.028 

± 

0.003 

5649.684 

1.1 

2.320 

± 

0.004 

2.899 

± 

0.004 

0.049 

± 

0.002 

5653.680 

0.9 

2.322 

± 

0.004 

2.890 

± 

0.004 

0.060 

± 

0.002 

(5664.715) 

(1.7) 

(2.331 

± 

0.005) 

(2.901 

± 

0.006) 

(-0.030 

± 

0.003) 

(5674.712) 

(1.7) 

(2.357 

± 

0.008) 

(2.838 

± 

0.011) 

(-0.127 

± 

0.004) 

(5684.642) 

(1.9) 

(2.332 

± 

0.005) 

(2.911 

± 

0.006) 

(0.011 

± 

0.003) 

5702.650 

1.5 

2.323 

± 

0.006 

2.884 

± 

0.007 

-0.067 

± 

0.003 

5712.657 

1.0 

2.330 

± 

0.005 

2.893 

± 

0.005 

0.012 

± 

0.003 

5864.959 

1.6 

2.280 

± 

0.004 

2.828 

± 

0.005 

0.026 

± 

0.002 

5888.932 

1.6 

2.286 

± 

0.004 

2.825 

± 

0.005 

0.017 

± 

0.002 

(5921.935) 

(1.7) 

(2.284 

± 

0.005) 

(2.854 

± 

0.006) 

(0.007 

± 

0.003) 

5924.982 

1.3 

2.287 

± 

0.005 

2.855 

± 

0.006 

0.009 

± 

0.003 

5947.799 

1.6 

2.288 

± 

0.004 

2.825 

± 

0.005 

0.026 

± 

0.002 

(5956.696) 

(2.0) 

(2.291 

± 

0.005) 

(2.838 

± 

0.006) 

(-0.016 

± 

0.003) 

(5975.687) 

(2.2) 

(2.411 

± 

0.006) 

(2.639 

± 

0.007) 

(-0.038 

± 

0.002) 

Notes.  HJD  is  the  Heliocentric  Julian  Day.  The  magnitudes  listed  in  the  QSO  A  and  B  columns  are  measured  relative  to 
the  comparison  stars.  The  magnitudes  in  the  (Stars)  column  are  the  mean  magnitudes  of  the  comparison  stars  for  that  epoch 
relative  to  their  mean  over  all  epochs.  The  light  curve  points  listed  in  parentheses  have  not  been  included  in  the  analysis. 
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In  addition,  here  we  present  data  from  the  g-band  monitoring 
program  of  SBS  0909  at  the  Liverpool  Telescope,  contempora¬ 
neous  with  the  r-band  data  set  published  in  Goicoechea  et  al. 
(2008)  and  thus  spanning  two  years  (2005  January  to  2007 
January).  We  used  the  RATCam  instrument  at  the  Liverpool 
Telescope  and  obtained  167  individual  exposures  (frames)  of 
100  or  200  s  each.  After  bias  subtraction,  overscan  region  trim¬ 
ming,  and  flat  fielding  of  the  images,  a  crowded-field  PSF  pho¬ 
tometry  pipeline  measures  instrumental  fluxes  for  bright  field 
stars  and  the  quasar  images.  We  then  transform  the  photometry 
to  the  SDSS  photometric  system,  correcting  the  instrumental 
fluxes  for  the  systematic  effects  of  color  and  inhomogeneity 
(see  Goicoechea  et  al.  2010).  The  transformation  pipeline  is 
only  applied  to  the  frames  in  which  the  signal-to-noise  ratio  of 
the  “c”  field  star,  measured  through  an  aperture  of  radius  equal 
to  twice  the  FWHM,  is  greater  than  100,  and  for  which  the  see¬ 
ing  is  less  than  2".  We  also  discard  frames  requiring  anomalous 
color  coefficients,  frames  which  produce  photometry  outliers, 
and  frames  in  which  the  quality  of  the  PSF  fits  to  the  quasar 
region  are  poor.  After  averaging  together  the  photometry  from 
individual  frames  obtained  on  the  same  night,  the  final  Liver¬ 
pool  Telescope  g-band  data  set  consists  of  43  epochs  of  SDSS 
magnitudes  with  average  uncertainties  of  0.016  mag  (image  A) 
and  0.036  mag  (image  B). 

The  new  r-band  photometry  from  the  Liverpool  Telescope  is 
listed  in  Table  2  and  the  complete  g-band  data  set  is  provided  in 
Table  3.  We  list  the  Liverpool  Telescope  light  curves  separately 
from  the  USNO  light  curves  because  the  two  data  sets  have 
different  photometric  calibrations:  the  Liverpool  data  frames 
have  been  calibrated  onto  an  absolute  system  using  the  absolute 
flux  of  a  reference  star,  while  the  USNO  measurements  are 
not  calibrated  to  a  standard  system.  By  presenting  the  two  data 
sets  separately,  we  preserve  the  original  photometric  system  of 
each  and  provide  transparency  of  origin  for  future  users  of  the 
data  sets. 

2.3.  Construction  of  Light  Curves  and  Difference  Light  Curves 

In  order  to  construct  light  curves  spanning  the  longest  possible 
time  baseline,  we  combined  all  the  r-band  USNO  and  Liverpool 
data  for  SBS  0909  along  with  one  R-band  epoch  obtained 
at  the  MDM  Observatory’s  Hiltner  2.4  m  telescope,  using 
the  1024  x  1024  “Templeton”  CCD  camera,  and  two  SDSS 
r-band  epochs  obtained  with  the  WIYN  Tip-Tilt  Module  at 
the  Wisconsin-Indiana-Yale-NOAO  (WIYN)  3.5  m  telescope. 
Since  these  epochs  were  not  contemporaneous  with  each  other 
or  our  USNO/Liverpool  light  curves,  we  were  unable  to  make 
an  empirical  measurement  of  any  magnitude  offsets  arising  from 
differences  between  the  detectors  and  filters.  We  accounted 
for  this  unknown  offset  by  applying  an  additional  0.02  mag 
of  uncertainty  to  the  time-delay  corrected  flux  ratio  for  these 
observations.  We  determined  the  magnitude  offset  between  our 
USNO  and  Liverpool  data  sets  (—14.108  ±  0.018)  by  making  a 
weighted  average  of  the  offset  A(USNO  —  LRT)  found  for  the 
seven  individual  nights  with  contemporaneous  observations. 

In  Figure  1  we  show  the  combined  r-band  light  curves  for 
images  A  and  B  from  all  of  our  data  sources.  The  dominant 
feature  in  the  light  curves  is  the  intrinsic  variability,  as  much 
as  0. 3-0.4  mag,  the  analysis  of  which  we  will  present  in 
a  future  paper.  Closer  inspection,  though,  of  the  time  range 
of  4500  <  HJD  -  2450000  <  5600  reveals  uncorrelated 
variability,  as  the  slope  of  the  increase  in  the  brightness  of 
image  B  over  this  time  period  is  notably  steeper  than  the  slope 
of  image  A’s  light  curve.  We  attribute  this  component  of  the 


Table  2 

SBS  0909+532  r  Light  Curves  from  Liverpool  Telescope 


HJD  -  2450000 
(days) 

QSO  A 
(mag) 

QSO  B 
(mag) 

4009.709 

16.379  ±0.010 

17.073  ±  0.013 

4011.697 

16.394  ±0.014 

17.019  ±0.018 

4015.729 

16.394  ±0.008 

17.053  ±0.010 

4017.705 

16.397  ±0.014 

17.042  ±0.018 

4028.658 

16.419  ±0.014 

17.040  ±0.018 

4029.668 

16.416  ±0.014 

17.156  ±0.018 

4030.708 

16.403  ±0.014 

17.117  ±0.018 

4031.654 

16.393  ±0.014 

17.097  ±0.018 

4049.600 

16.436  ±0.014 

17.031  ±0.018 

4057.601 

16.438  ±0.014 

17.031  ±0.018 

4062.564 

16.405  ±0.014 

17.078  ±0.018 

4063.566 

16.394  ±0.014 

17.094  ±0.018 

4069.536 

16.395  ±0.014 

17.106  ±0.018 

4070.570 

16.453  ±0.014 

17.043  ±0.018 

4074.556 

16.437  ±0.014 

17.062  ±0.018 

4077.576 

16.489  ±0.014 

17.021  ±0.018 

4084.533 

16.415  ±0.014 

17.109  ±0.018 

4085.553 

16.411  ±0.014 

17.132  ±0.018 

4091.483 

16.381  ±0.014 

17.131  ±0.018 

4101.453 

16.408  ±0.014 

17.159  ±0.018 

4105.604 

16.394  ±0.014 

17.160  ±  0.018 

4109.528 

16.419  ±0.014 

17.166  ±0.018 

4110.523 

16.410  ±0.014 

17.229  ±0.018 

4111.431 

16.412  ±0.014 

17.121  ±0.018 

4115.450 

16.440  ±0.014 

17.1 10  ±  0.018 

4117.563 

16.412  ±0.014 

17.198  ±0.018 

4120.409 

16.370  ±0.014 

17.207  ±0.018 

4122.508 

16.414  ±0.014 

17.187  ±0.018 

4123.439 

16.438  ±0.014 

17.135  ±0.018 

4124.600 

16.430  ±0.014 

17.185  ±0.018 

5474.697 

16.458  ±0.014 

16.992  ±0.018 

5479.672 

16.469  ±0.014 

17.015  ±0.018 

5486.683 

16.465  ±0.014 

16.966  ±0.018 

5488.734 

16.437  ±0.014 

17.006  ±0.018 

5492.716 

16.445  ±0.014 

17.008  ±0.018 

5494.700 

16.478  ±0.014 

16.968  ±0.018 

5496.656 

16.468  ±0.014 

17.002  ±0.018 

5500.743 

16.429  ±0.014 

17.043  ±0.018 

5501.747 

16.446  ±0.014 

17.023  ±0.018 

5506.697 

16.439  ±0.014 

17.010±  0.018 

5507.611 

16.431  ±0.014 

17.025  ±0.018 

5510.653 

16.442  ±0.014 

17.009  ±0.018 

5511.604 

16.453  ±0.014 

16.989  ±0.018 

5512.613 

16.457  ±0.014 

17.005  ±0.018 

5513.700 

16.449  ±0.014 

17.002  ±0.018 

5514.594 

16.466  ±0.014 

16.967  ±0.018 

5516.698 

16.454  ±0.014 

16.992  ±0.018 

5517.599 

16.449  ±0.014 

17.000  ±0.018 

5518.584 

16.458  ±0.014 

16.995  ±0.018 

5519.664 

16.442  ±0.014 

16.995  ±0.018 

5520.677 

16.449  ±0.014 

16.987  ±0.018 

5522.696 

16.453  ±0.014 

17.004  ±0.018 

5524.623 

16.446  ±0.014 

17.003  ±0.018 

5531.640 

16.470  ±0.014 

16.954  ±0.018 

5553.519 

16.471  ±0.014 

16.957  ±0.018 

5554.527 

16.452  ±0.014 

16.951  ±0.018 

5562.640 

16.456  ±0.014 

16.974  ±0.018 

5563.472 

16.462  ±0.014 

16.947  ±0.018 

5564.441 

16.418  ±0.014 

17.010  ±0.018 

5566.511 

16.437  ±0.014 

17.003  ±0.018 

5567.576 

16.435  ±0.014 

17.003  ±0.018 

5568.499 

16.437  ±0.014 

16.980  ±0.018 

5569.630 

16.437  ±0.014 

16.979  ±0.018 

5571.418 

16.423  ±0.014 

16.975  ±0.018 

5572.424 

16.419  ±0.014 

17.000  ±0.018 
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Table  2 

(Continued) 


HJD  -  2450000 
(days) 

QSO  A 
(mag) 

QSO  B 
(mag) 

5573.510 

16.437 

± 

0.014 

16.987 

±0.018 

5579.507 

16.415 

± 

0.014 

17.007 

±0.018 

5580.429 

16.438 

± 

0.014 

16.988 

±0.018 

5581.494 

16.430 

± 

0.014 

17.005 

±0.018 

5599.429 

16.429 

± 

0.014 

17.024 

±0.018 

5600.475 

16.424 

± 

0.014 

17.013 

±0.018 

5605.447 

16.408 

± 

0.014 

17.025 

±0.018 

5607.438 

16.414 

± 

0.014 

17.018 

±0.018 

5609.423 

16.418 

± 

0.014 

17.051 

±0.018 

5610.400 

16.431 

± 

0.014 

17.002 

±0.018 

5611.364 

16.430 

± 

0.014 

17.015 

±0.018 

5612.408 

16.408 

± 

0.014 

17.032 

±0.018 

5615.365 

16.427 

± 

0.014 

17.002 

±0.018 

5616.452 

16.437 

± 

0.014 

17.013 

±0.018 

5617.366 

16.430 

± 

0.014 

17.006 

±0.018 

5619.365 

16.432 

± 

0.014 

17.011 

±0.018 

5620.384 

16.437 

± 

0.014 

17.007 

±0.018 

5621.368 

16.421 

± 

0.014 

17.014 

±0.018 

5622.368 

16.424 

± 

0.014 

17.012 

±0.018 

5623.360 

16.437 

± 

0.014 

17.018 

±0.018 

5644.407 

16.434 

± 

0.014 

17.018 

±0.018 

5645.372 

16.445 

± 

0.014 

17.029 

±0.018 

5646.369 

16.421 

± 

0.014 

17.017 

±0.018 

5647.374 

16.432 

± 

0.014 

17.028 

±0.018 

5649.377 

16.435 

± 

0.014 

17.011 

±0.018 

5650.386 

16.423 

± 

0.014 

17.018 

±0.018 

5651.375 

16.436 

± 

0.014 

17.000 

±0.018 

5652.398 

16.431 

± 

0.014 

17.022 

±0.018 

5653.377 

16.430 

± 

0.014 

17.003 

±0.018 

5655.381 

16.428 

± 

0.014 

17.004 

±0.018 

5656.465 

16.419 

± 

0.014 

17.047 

±0.018 

5657.373 

16.412 

± 

0.014 

17.022 

±0.018 

5658.403 

16.417 

± 

0.014 

17.029 

±0.018 

5660.426 

16.430 

± 

0.014 

17.031 

±0.018 

5661.373 

16.437 

± 

0.014 

17.023 

±0.018 

5662.375 

16.428 

± 

0.014 

17.015 

±0.018 

5663.389 

16.409 

± 

0.014 

17.042 

±0.018 

5664.379 

16.426 

± 

0.014 

17.029 

±0.018 

5689.395 

16.424 

± 

0.014 

17.031 

±0.018 

5690.390 

16.416 

± 

0.014 

17.063 

±0.018 

5691.437 

16.422 

± 

0.014 

17.024 

±0.018 

5692.454 

16.431 

± 

0.014 

17.015 

±0.018 

5693.397 

16.418 

± 

0.014 

17.023 

±0.018 

5694.416 

16.424 

± 

0.014 

17.026 

±0.018 

5696.411 

16.419 

± 

0.014 

16.998 

±0.018 

5698.455 

16.416 

± 

0.014 

17.030 

±0.018 

5699.437 

16.435 

± 

0.014 

17.001 

±0.018 

5701.412 

16.432 

± 

0.014 

16.992 

±0.018 

5703.432 

16.417 

± 

0.014 

17.024 

±0.018 

5705.403 

16.414 

± 

0.014 

17.014 

±0.018 

5706.490 

16.415 

± 

0.014 

17.008 

±0.018 

5710.406 

16.430 

± 

0.014 

17.009 

±0.018 

5712.403 

16.431 

± 

0.014 

16.997 

±0.018 

5713.451 

16.423 

± 

0.014 

17.019 

±0.018 

5714.403 

16.413 

± 

0.014 

17.001 

±0.018 

Notes.  HJD  is  the  Heliocentric  Julian  Day.  The  magnitudes  listed 
have  been  calibrated  to  the  SDSS  photometric  system  using  the  flux 
of  the  “b”  star  in  the  field  of  SBS  0909  (see  Kochanek  et  al.  1997). 

variability  to  microlensing  by  stars  in  the  lens  galaxy.  A  careful 
examination  of  the  lightcurves  and  Tables  1  and  2  will  reveal 
that  our  data  cadence  is  somewhat  variable.  In  the  first  few 

seasons  there  are  many  periods  with  several  observations  per 


Table  3 

SBS  0909+532  g  Light  Curves  from  Liverpool  Telescope 


HJD  -  2450000 
(days) 

QSO  A 
(mag) 

QSO  B 
(mag) 

3649.708 

17.002 

± 

0.017 

18.053 

± 

0.038 

3656.720 

17.008 

± 

0.017 

18.001 

± 

0.038 

3676.666 

17.025 

± 

0.010 

17.998 

± 

0.022 

3677.665 

17.011 

± 

0.012 

17.997 

± 

0.027 

3684.674 

16.980 

± 

0.010 

18.079 

± 

0.022 

3700.678 

17.005 

± 

0.010 

18.028 

± 

0.022 

3709.702 

17.001 

± 

0.012 

18.063 

± 

0.027 

3711.699 

17.001 

± 

0.017 

18.058 

± 

0.038 

3714.595 

17.010 

± 

0.010 

18.054 

± 

0.022 

3729.504 

17.057 

± 

0.017 

17.949 

± 

0.038 

3731.561 

17.020 

± 

0.008 

18.006 

± 

0.019 

3784.540 

17.001 

± 

0.017 

18.017 

± 

0.038 

3786.463 

16.992 

± 

0.017 

18.044 

± 

0.038 

3787.458 

16.989 

± 

0.017 

18.024 

± 

0.038 

3789.450 

16.986 

± 

0.017 

18.024 

± 

0.038 

3790.423 

16.990 

± 

0.017 

18.003 

± 

0.038 

3791.424 

16.998 

± 

0.017 

18.017 

± 

0.038 

3793.411 

16.989 

± 

0.017 

18.011 

± 

0.038 

3798.447 

16.994 

± 

0.017 

17.951 

± 

0.038 

3802.407 

16.984 

± 

0.017 

17.985 

± 

0.038 

3813.396 

16.980 

± 

0.017 

17.975 

± 

0.038 

3814.426 

16.968 

± 

0.017 

17.992 

± 

0.038 

3815.386 

16.973 

± 

0.017 

17.964 

± 

0.038 

3817.440 

16.953 

± 

0.017 

18.033 

± 

0.038 

3818.427 

16.958 

± 

0.017 

18.016 

± 

0.038 

3827.521 

16.958 

± 

0.017 

18.051 

± 

0.038 

3830.475 

16.971 

± 

0.017 

18.021 

± 

0.038 

3845.429 

16.955 

± 

0.017 

18.020 

± 

0.038 

3854.406 

16.966 

± 

0.017 

18.029 

± 

0.038 

3855.491 

16.951 

± 

0.017 

18.053 

± 

0.038 

3864.456 

16.975 

± 

0.017 

18.043 

± 

0.038 

3875.420 

16.967 

± 

0.017 

18.066 

± 

0.038 

3887.418 

16.958 

± 

0.017 

18.037 

± 

0.038 

3888.420 

16.934 

± 

0.017 

18.101 

± 

0.038 

3890.421 

16.947 

± 

0.017 

18.065 

± 

0.038 

4009.714 

16.996 

± 

0.017 

18.043 

± 

0.038 

4030.710 

16.991 

± 

0.017 

18.100 

± 

0.038 

4031.656 

17.036 

± 

0.017 

17.989 

± 

0.038 

4032.647 

17.027 

± 

0.017 

18.029 

± 

0.038 

4062.566 

17.027 

± 

0.017 

18.104 

± 

0.038 

4085.555 

17.031 

± 

0.017 

18.131 

± 

0.038 

4122.509 

17.018 

± 

0.017 

18.205 

± 

0.038 

4124.602 

17.008 

± 

0.017 

18.240 

± 

0.038 

Notes.  HJD  is  the  Heliocentric  Julian  Day.  The  magnitudes  listed 
have  been  calibrated  to  the  SDSS  photometric  system  using  the  flux 
of  the  “b”  star  in  the  field  of  SBS  0909  (see  Kochanek  et  al.  1997). 

week,  but  much  of  the  USNO  data  is  at  a  cadence  of  one  to  two 
observations  per  month.  The  strength  of  our  analysis  technique 
is  that  while  the  sparsely  sampled  periods  do  not  constrain 
the  time  delay  directly,  their  long  time  baseline  provides  very 
strong  constraints  on  the  microlensing  model,  thereby  severely 
limiting  the  number  of  trial  light  curves  from  our  Monte  Carlo 
simulation  with  adequate  fits  to  the  data  during  the  densely 
sampled  intervals  (see  Section  3.1). 

We  display  the  shorter  observed  g-band  light  curves  for 

SBS  0909,  which  are  composed  entirely  of  Liverpool  Telescope 
data,  in  Figure  2.  The  intrinsic  variability  in  the  g-band  is  similar 
to  that  in  the  /'-band  during  the  same  time  period.  The  g-band 
curves  exhibit  considerably  more  scatter  than  the  r-band  curves, 
particularly  for  image  B .  The  increased  scatter  is  due  primarily  to 
the  lower  quasar  flux  in  the  g  band  and  the  increased  difficulty  in 
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Figure  1.  Composite  r  band  light  curves  for  SBS  0909  images  A  (top)  and  B 
(bottom)  including  measurements  from  MDM  Observatory  (stars),  the  WIYN 
3.5  m  telescope  (diagonal  crosses),  the  Liverpool  Telescope  (triangles),  and 
USNO  (squares).  The  measurements  for  image  A  have  been  offset  by  +0.2  mag 
to  minimize  empty  space  in  the  plot  area.  The  light  curve  of  image  B  exhibits 
a  substantially  steeper  slope  over  the  time  period  4500  <  HJD  —  2450000  < 
5200  days. 


Figure  2.  Observed  g  band  light  curves  for  SBS  0909  images  A  (top  panel)  and 
B  (bottom  panel)  from  the  Liverpool  Telescope.  The  g-band  light  curves  exhibit 
similar  intrinsic  variability  to  the  r  band  light  curves  over  the  same  period  of 
time,  although  with  increased  scatter  due  to  the  lower  quasar  flux  in  g  and  poorer 
observing  conditions  on  some  occasions. 


extracting  bluer  g-band  fluxes  under  poor  observing  conditions. 
Although  the  increased  scatter  impedes  visual  identification 
of  microlensing  variability  in  the  g  band  light  curves,  the 
measurements  are  still  valuable  as  a  complementary  data  set 
for  a  simultaneous  two-band  microlensing  analysis. 

The  microlensing  variability  in  the  g  and  r  band  light  curves 
is  displayed  in  Figure  3.  These  time-delay-shifted  difference 
light  curves  are  generated  by  shifting  the  light  curve  of  the  less 
variable  image  (image  A,  in  both  bands)  by  the  system’s  50  day 
time  delay  A tAB  and  then  performing  a  linear  interpolation  of 
image  A’s  shifted  light  curve  to  generate  a  set  of  photometric 
measurements  at  the  same  epochs  of  observation  as  those  in 
image  B ’s  (unshifted)  light  curve.  We  discard  any  data  points  that 


HJD  -  2450000  (days) 


HJD  -  2450000  (days) 


Figure  3.  Difference  light  curves  in  the  g  (top  panel)  and  r  bands  (bottom  panel) 
for  SBS  0909,  shown  with  an  example  of  a  simulated  light  curve  from  our 
Monte  Carlo  simulations  that  is  a  good  fit  to  the  observations.  To  construct  the 
light  curves,  image  A’s  data  has  been  shifted  by  A  t\B  =  tA  —  ts  =  5  0  days. 
Significant  uncorrelated  variability  is  apparent  in  the  r  band. 

were  interpolated  in  the  interseason  gaps.  Finally,  we  subtract 
the  light  curve  of  image  B  from  the  shifted  and  interpolated 
light  curve  of  image  A.  We  do  not  apply  any  corrections  for  lens 
galaxy  reddening  to  the  light  curves  or  difference  light  curve. 
Rather,  we  allow  for  a  0.5  mag  systematic  uncertainty  in  the 
intrinsic  value  flux  ratios  in  our  simulations,  which  accounts  for 
the  uncertainty  in  the  intrinsic  flux  ratio  from  both  microlensing 
and  dust  extinction  in  the  lens  galaxy  (Mediavilla  et  al.  2011). 
The  resulting  difference  light  curves  are  shown  in  Figure  3  for 
both  g  and  r.  As  expected,  the  microlensing  signal  is  difficult 
to  identify  in  the  shorter  and  noisier  g-band  difference  light 
curve.  However,  the  uncorrelated  variability  in  the  r  band  is 
quite  clear.  We  observe  a  steep  fall  and  rise  in  A m  A  —  Am  B  in 
the  time  interval  4000  <  HJD  -  2450000  <  5200,  followed 
by  a  second,  smaller-amplitude  oscillation  in  the  time  interval 
5500  <  HJD  -  2450000  <  5900. 

3.  MONTE  CARLO  MICROLENSING 
AND  TIME  DELAY  ANALYSIS 

As  discussed  in  Section  1,  we  wish  to  use  our  combined 
r  and  g  band  light  curves  to  derive  a  more  precise  time 
delay  for  SBS  0909  as  well  as  to  place  new  constraints  on 
the  quasar  continuum  source  size  and  structure.  We  do  so 
here  utilizing  the  alternative  implementation  of  the  Kochanek 
(2004)  microlensing  analysis  techniques  presented  in  Morgan 
et  al.  (2008,  2012)  in  which  the  time  delay  and  source  size 
of  a  lensed  quasar  are  simultaneously  determined  through  a 
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Bayesian  analysis  of  Monte  Carlo  microlensing  simulations  for 
a  range  of  trial  time  delays.  Our  approach  here  differs  only  in 
that  we  analyze  two  optical  light  curves  (r  and  g  band)  for  each 
quasar  image  rather  than  an  optical  light  curve  and  an  X-ray 
light  curve  for  each  image.  We  carry  out  the  microlensing  and 
time  delay  analysis  in  two  major  steps:  first,  we  analyze  Monte 
Carlo  simulations  of  the  longer,  better-sampled  light  curves  (in 
this  case,  the  r  band)  for  a  range  of  time  delays  and  source  sizes. 
Then,  we  simultaneously  analyze  Monte  Carlo  simulations  of 
the  g  and  r  light  curves  incorporating  the  best  time  delay  value 
and  the  continuum  source  size  distribution  produced  in  the  first 
stage.  The  final  results  for  the  r-  and  g-band  continuum  source 
sizes  are  derived  in  this  second  step. 

3.1.  Time  Delay  Analysis  of  r -band  Light  Curves 

We  analyze  only  the  r  band  light  curves  for  the  time  delay 
analysis,  as  the  r-band  monitoring  spans  a  longer  length  of  time 
and  includes  the  largest-amplitude  flux  variation.  To  begin  the 
analysis,  we  generate  r  band  light  curve  pairs  in  which  image 
A’s  (the  less  variable  image)  light  curve  is  shifted  by  a  set  of  trial 
time  delays,  AtA  B  =  tA  —  tB,  according  to  the  procedure  outlined 
in  Morgan  et  al.  (2008).  Our  set  of  trial  time  delays  spans  the 
range  —70  days  -<)  A tAB  f  70  days  in  time  steps  of  one  day.  The 
shifting  of  the  light  curves  by  the  trial  time  delays  must  take  into 
account  two  issues.  First,  as  explained  in  Morgan  et  al.  (2008, 
2012),  each  A/B  light  curve  pair  must  have  measurements  for 
the  same  dates,  necessitating  interpolation  of  the  shifted  image 
A  measurements  when  they  fall  in  the  middle  of  an  observing 
season  and  extrapolation  of  image  A  measurements  when  the 
shifted  data  points  fall  in  interseason  gaps.  We  used  linear 
interpolation  for  epochs  falling  in  the  middle  of  an  observing 
season  and  permitted  10  days  of  linear  extrapolation  for  points 
in  interseason  gaps,  increasing  the  uncertainties  for  successive 
extrapolated  points  as  described  in  Morgan  et  al.  (2012).  The 
second  issue  is  that  all  trial  light  curves  must  have  the  same 
number  of  epochs,  requiring  the  truncation  of  light  curves  for 
trial  time  delays  shorter  than  60  days.  After  all  interpolations, 
extrapolations,  and  truncations,  the  shifted  light  curves  used  for 
our  analysis  each  contained  158  epochs. 

Next,  for  each  trial  time  delay,  we  use  the  Monte  Carlo  method 
of  Kochanek  (2004;  also  see  Poindexter  &  Kochanek  2010)  to 
generate  large  numbers  of  light  curves  from  microlensing  mag¬ 
nification  patterns  for  random  combinations  of  effective  velocity 
ve  between  quasar,  lens  galaxy,  and  observer,  mean  microlens 
mass,  and  lens  galaxy  stellar  mass/dark  matter  fraction.  We 
then  fit  the  simulated  light  curves  to  the  observed  r  band  light 
curve  shifted  by  the  trial  time  delay.  To  construct  the  magnifi¬ 
cation  patterns,  we  must  utilize  physical  models  of  the  macro¬ 
scopic  (strong)  lensing,  microlens  mass  function,  and  accretion 
disk  surface  brightness  profile.  We  used  the  lensmodel  software 
package  (Keeton  2001)  to  generate  a  sequence  of  strong  lensing 
models  for  the  SBS  0909  system  over  a  range  of  mass  contri¬ 
butions  from  the  dark  matter  and  stellar  components  of  the  lens 
galaxy.  Each  model  is  a  sum  of  a  de  Vaucouleurs  component, 
representing  the  stellar  content  of  the  lens  galaxy,  and  a  Navarro 
et  al.  (1996,  NFW)  component,  representing  the  galaxy’s  dark 
matter  halo  and  concentric  with  the  de  Vaucouleurs  component. 
We  ran  our  first  realization  of  the  Monte  Carlo  simulation  using 
a  model  sequence  whose  coefficients  best  reproduce  the  lens 
galaxy  model  of  Lehar  et  al.  (2000),  but  we  eventually  used  a 
range  of  model  sequences  at  both  the  Lehar  and  Sluse  et  al. 
(2012)  positions.  A  model  sequence  contains  10  models  span¬ 
ning  0.1  ^  /m/l  ^  1.0  in  steps  of  0.1,  where  /m/l  represents 


the  ratio  of  the  mass  of  the  stellar  component  to  its  mass  in  a  uni¬ 
form  mass-to-light  ratio  model.  For  the  stellar  (microlens)  mass 
function,  we  use  a  power  law,  dN/dM  oc  M~13,  with  a  ratio 
of  maximum-to-minimum  mass  ratio  of  50;  this  function  rea¬ 
sonably  approximates  the  Galactic  disk  mass  function  of  Gould 
(2000).  We  model  the  quasar’s  accretion  disk  as  a  face-on,  thin 
disk  radiating  as  a  blackbody  with  a  power-law  temperature  pro¬ 
file  T  oc  R  7/4 .  Our  model  matches  the  outer  regions  of  the  thin 
disk  model  of  Shakura  &  Sunyaev  (1973),  but  we  neglect  the 
drop  in  temperature  in  the  center  due  to  the  inner  edge  of  the  disk 
and  the  correction  factor  from  general  relativity  to  avoid  intro¬ 
ducing  additional  parameters.  Provided  the  disk  sizes  we  obtain 
are  significantly  larger  than  the  radius  of  the  inner  disk  edge, 
these  simplifications  introduce  insignificant  uncertainties  rela¬ 
tive  to  those  associated  with  other  parameters  (Dai  et  al.  2010). 
With  these  parameters  we  create  40  independent  magnification 
patterns  for  each  quasar  image  for  each  of  the  10  different  strong 
lens  models,  using  the  method  described  in  the  Appendix  of 
Kochanek  (2004).  The  patterns  are  8192  x  8192  images  repre¬ 
senting  20 (Re)  x  20 (Re),  where  (  Re)  is  the  Einstein  radius  for 
the  mean  microlens  mass  ( M )  projected  into  the  source  plane, 
yielding  a  pixel  scale  of  1.1  x  1014((M) /M0)1//2  cm.  The  outer 
dimensions  and  pixel  scale  are  chosen  to  be  sufficiently  large  to 
representatively  sample  the  magnification  distribution  and  suf¬ 
ficiently  small  to  adequately  resolve  the  accretion  disk  in  the 
g-band  simulations  (see  Section  3.2). 

For  the  time  delay  analysis  we  carry  out  107  realizations  of 
the  r  band  light  curve  from  each  of  the  400  sets  of  magnification 
patterns  and  for  each  of  the  141  trial  time  delays.  We  randomly 
select  an  initial  position  and  effective  velocity  for  the  source 
trajectory  from  their  prior  distributions  under  the  assumption 
that  these  variables  are  independent  and  uniformly  distributed. 
For  computational  simplicity,  we  neglect  the  motion  of  the  stars 
within  the  lens  galaxy  and  describe  the  observer’s  motion  as  the 
projection  of  the  cosmic  microwave  background  dipole  velocity 
onto  the  lens  plane,  as  done  by  Kochanek  (2004).  We  compare 
the  simulated  light  curves  to  the  observed  light  curves  and 
calculate  the  goodness-of-fit  (/2)  statistics  for  each,  discarding 
trials  with  a  /2  statistic  per  degree  of  freedom  (x2/v)  greater 
than  3. 

As  we  discussed  in  Section  1 ,  the  time  delay  measurement  in 
SBS  0909  is  made  more  challenging  by  the  uncertainties  in  the 
lens  galaxy  photometric  model.  Since  Lehar  et  al.  (2000)  and 
Sluse  et  al.  (2012)  disagree  about  the  astrometric  position  and 
effective  radius  0'eff)  of  the  lens  galaxy,  we  took  steps  to  ensure 
the  robustness  of  our  time  delay  measurement  in  the  resulting 
degenerate  model  parameter  space.  First,  we  created  a  series  of 
additional  model  sequences  for  a  range  of  lens  galaxy  effective 
radii  (reff)  at  the  astrometric  positions  as  determined  by  both 
Sluse  et  al.  (2012)  and  Lehar  et  al.  (2000).  Each  of  these  model 
sequences  employs  a  two  component  (de  Vaucouleurs— NFW) 
lens  galaxy,  where  we  vary  the  stellar  mass  fraction  as  described 
above,  imposing  ellipticity  constraints  on  the  lens  galaxy  from 
Sluse  et  al.  (2012)  in  the  former  case,  and  a  small  ellipticity  of 
1  —  b/a  —  0  ±  0.08  in  the  latter  case.  Then,  for  each  new  model 
sequence,  we  generated  a  new  set  of  magnification  patterns 
with  which  we  repeated  the  full-scale  Monte  Carlo  time  delay 
analysis  described  above.  For  completeness,  we  also  repeated 
this  procedure  for  a  Singular  Isothermal  Ellipsoid  model  at  both 
the  Sluse  et  al.  (2012)  and  the  Lehar  et  al.  (2000)  positions. 

Following  the  completion  of  the  Monte  Carlo  light  curve 
simulations  for  all  model  sequences,  we  aggregated  the  results 
from  all  trials  and  performed  a  Bayesian  analysis  of  the  y1 
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Figure  4.  Posterior  probability  distribution  for  the  time  delay  in  SBS  0909.  The 
portion  of  the  distribution  for  delays  —70  days  <  A Iab  <  Odays  is  not  shown 
because  the  probability  in  that  section  is  essentially  zero.  Our  result  for  the  time 
delay,  A tAs  =  50^  days,  where  B  leads  A,  agrees  with  the  previous  result  from 
Goicoechea  et  al.  (2008),  but  is  more  precise. 

statistics  of  the  light  curve  fits.  We  formally  marginalize  over 
the  microlensing  variables  of  source  size,  microlens  mass, 
and  stellar  mass  fraction  used  to  construct  the  magnification 
patterns,  as  well  as  the  effective  velocities  from  the  Monte  Carlo 
simulations,  to  calculate  the  posterior  probability  density  for  the 
time  delay  of  SBS  0909.  By  aggregating  the  results  from  the 
model  sequences  at  both  the  Lehar  et  al.  (2000)  and  Sluse  et  al. 
(2012)  positions,  we  effectively  marginalize  over  uncertainties 
in  the  lens  galaxy  photometric  fits  as  well.  The  trials  with  high 
X1  which  were  thrown  out  would  not  contribute  significantly  to 
the  Bayesian  integrals,  so  removing  them  does  not  affect  our 
results. 

We  present  the  posterior  probability  density  for  the  time 
delay  At  a B  of  SBS  0909  resulting  from  our  Bayesian  analysis 
of  the  full  set  of  model  sequences  in  Figure  4.  For  ease  of 
viewing,  we  show  only  the  portion  of  the  distribution  for 
Odays  <  A  1AB  <  70 days;  the  values  of  the  probability  density 
for  negative  time  delays  were  essentially  zero.  The  time  delay 
distribution  is  narrowly  peaked,  with  a  median  A tAs  —  50  days 
(so  image  B  leads  image  A)  and  a  68%  confidence  interval  of 
46  days  <  AtA  <  52  days.  Our  result  is  consistent  with,  but 
more  precise  than,  the  previous  time  delay  measurement  A tAB  = 
49  ±  6  days  by  Goicoechea  et  al.  (2008),  based  on  the  first  two 
seasons  of  r-band  monitoring  data  used  in  this  study.  We  use 
this  revised  time  delay  for  the  remainder  of  our  analysis,  while 
acknowledging  that  we  were  unable  to  estimate  the  influence  of 
all  possible  systematic  errors  in  our  measurement  technique.  So, 
as  is  the  case  with  any  sophisticated  measurement,  it  is  possible 
that  we  underestimated  the  influence  of  systematic  errors,  and 
these  unknown  systematics  may  contribute  to  the  discrepancy 
we  describe  in  the  next  paragraph  and  explore  in  Section  4. 

In  Figure  5,  we  illustrate  the  influence  of  the  discrepancies  in 
the  lens  galaxy  photometric  fits  from  the  literature.  We  display 
contour  plots  of  the  time  delay  predicted  by  the  SBS  0909 
lens  models  as  a  function  of  the  effective  radius  req-  and 


fractional  mass  of  the  de  Vaucouleurs  component  relative  to 
constant  M/L  model,  /m/l,  for  lens  galaxy  astrometry  from 
both  Sluse  et  al.  (2012)  and  Lehar  et  al.  (2000).  As  expected, 
the  predicted  time  delay  is  longer  for  the  more  compact  lens 
models  with  larger  stellar  mass  components,  but  we  also  see 
that  the  delays  are  systematically  longer  for  a  galaxy  located  at 
the  Lehar  et  al.  (2000)  position  than  at  the  Sluse  et  al.  (2012) 
position.  This  difference  arises  because  the  lens  galaxy  in  the 
Lehar  et  al.  (2000)  fit  is  closer  to  image  A  than  in  the  Sluse 
et  al.  (2012)  fit,  yielding  a  larger  gravitational  delay  in  image 
A.  In  any  case,  neither  the  fiducial  Lehar  et  al.  (2000)  model 
nor  the  fiducial  Sluse  et  al.  (2012)  model  can  reproduce  our  new 
measurement  of  the  time  delay  AtAB  =  50+24  days,  but  galaxy 
models  at  the  Lehar  position  and  on  the  small  end  of  the  band  of 
uncertainty  in  the  Lehar  effective  radius  measurement  (reff  = 
1"58  ±  O'.' 9)  yield  delays  that  are  easily  consistent  with  our 
new  time  delay  measurement.  On  the  other  hand,  reproducing 
our  new  delay  measurement  with  models  at  the  Sluse  et  al. 
(2012)  position  requires  a  lens  galaxy  that  is  significantly 
smaller  than  the  already  much  more  compact  galaxy  in  Sluse 
measurement  (reff  =  0'.'54  ±  0'.'02).  Fortunately,  our  Monte 
Carlo  microlensing  simulation  is  sufficiently  realistic  as  to  be 
sensitive  to  the  differences  between  intrinsic  and  microlensing 
variability,  despite  the  uncertainties  in  the  macroscopic  lens 
model,  since  for  each  model  sequence  we  sample  a  wide  range 
of  the  stellar-to-total  convergence  ratio  k^/k.  with  significant 
overlap  between  the  macro  model  sequences. 

3.2.  Simultaneous  Dual-band  Microlensing  Analysis 

Our  dual-band  ( r  and  g)  microlensing  analysis  follows  the 
method  of  the  simultaneous  optical  and  X-ray  analyses  by  Dai 
et  al.  (2010)  and  Morgan  et  al.  (2012).  We  first  shift  the  r  and 
g  band  light  curves  by  the  new  time  delay,  A tAB  =  50  days,  in  the 
same  manner  as  we  employed  to  construct  the  difference  light 
curves  in  Section  2.3.  Using  the  same  magnification  patterns 
from  the  time  delay  analysis,  we  carry  out  107  simulations  of 
the  /--band  light  curve  for  each  of  the  sets  of  magnification 
patterns  and  discard  solutions  for  which  X2/v  >  2.5.  We  saved 
all  the  physical  parameters  from  the  surviving  r  band  light  curve 
fits.  We  then  attempted  to  fit  the  g-band  light  curve  using  the 
trajectories  from  the  best  r  band  fits  for  a  new  grid  of  source 
sizes,  and  we  compute  the  joint  x2/v  f°r  the  combined  r  and 
g  band  fits.  A  sample  simulated  r-  and  g-band  difference  light 
curve  which  is  a  best  fit  to  the  observed  data  is  shown  in  Figure  3 . 

We  calculate  posterior  probability  density  distributions  for 
the  accretion  disk  sizes  in  the  g-  and  /--hands  and  the  lens  galaxy 
stellar  mass  fraction  (/m/l)  by  performing  a  Bayesian  analysis 
on  the  combined  set  of  r-  and  g-band  solutions.  We  note  that 
our  simulations  are  carried  out  in  Einstein  units,  where  source 
sizes  and  velocity  are  scaled  by  ({M/Mq})1^2  and  denoted 
by  rs  and  ve ,  respectively.  To  obtain  the  probability  density 
for  the  true,  unsealed  physical  source  size  P(rs)  from  that 
for  the  scaled  source  size,  P(rs),  we  combine  P(rs)  with  the 
probability  density  for  the  scaled  effective  velocity,  P(ve),  and 
a  statistical  model  (i.e.,  a  prior)  for  the  true  effective  source 
velocity,  P(ve),  in  our  analysis.  We  construct  P(v,.)  using  the 
method  described  in  Kochanek  (2004).  For  that  purpose  we  use 
the  peculiar  velocity  estimates  for  the  redshifts  of  SBS  0909 
and  the  lens  galaxy  from  the  models  presented  in  Mosquera 
&  Kochanek  (2011)  and  estimate  the  velocity  dispersion  of 
the  lens  galaxy  from  its  Einstein  radius,  assuming  the  galaxy 
is  a  singular  isothermal  sphere  with  relaxed  dynamics,  which 
Treu  &  Koopmans  (2004)  and  Bolton  et  al.  (2008)  show  is  a 


8 


The  Astrophysical  Journal,  774:69  (llpp),  2013  September  1 


Hainline  et  al. 


Figure  5.  Contours  of  time  delay  A t&B  predicted  by  a  grid  of  de  Vaucouleurs  +  NFW  lens  models  as  a  function  of  the  de  Vaucouleurs  fraction  /m/l  and  effective 
radius  for  the  astrometric  solutions  of  Lehar  et  al.  (2000;  left  panel)  and  Sluse  et  al.  (2012;  right  panel).  The  contour  spacing  is  5  days.  For  each  effective  radius, 
the  base  model  is  a  de  Vaucouleurs-only  model  (/m/l  =  1-0)  in  which  the  mass  parameter  is  allowed  to  vary  to  yield  the  best  fit.  At  each  the  /m/l  parameter 
(see  Section  3.1)  varies  in  uniform  steps  between  1.0  and  0.1  relative  to  the  base  model. 


good  approximation.  As  a  final  step  we  must  correct  the  scale 
radius  for  the  disk’s  inclination  i  by  multiplying  by  (cos  /)'l/2, 
which  is  necessary  because  we  have  assumed  a  face-on  disk 
in  our  simulations  and  microlensing  amplitudes  depend  on  the 
projected  area  of  a  source  rather  than  the  shape. 

4.  RESULTS  AND  DISCUSSION 

In  the  top  panel  of  Figure  6,  we  show  the  posterior  probability 
distribution  for  the  physical  size  of  the  quasar's  accretion  disk 
in  the  observed-frame  r  band  resulting  from  our  two-band 
Bayesian  microlensing  analysis.  Hereafter,  we  state  all  sizes 
in  terms  of  the  thin  disk  scale  radius,  rs,  defined  as  the  radius  at 
which  the  disk  temperature  matches  the  rest-frame  wavelength 
of  the  filter  used  in  our  monitoring  observations,  kT  =  /;c//,rest 
(for  r- ha nd  monitoring  of  SBS  0909,  Arest  =  2620  A;  for 
g-band,  /,resj  =  2020  A).  This  can  be  converted  to  a  half- 
light  radius  using  the  relation  r\/2  =  2AArs.  We  note  that  the 
distribution  shown  and  all  numerical  quantities  in  the  discussion 
which  follows  have  been  corrected  for  disk  inclination  i  through 
multiplication  by  a  factor  of  (cos/)'1/2,  assuming  /  =  60° 
(corresponding  to  the  expectation  value  of  a  random  distribution 
of  disk  inclinations).  The  median  of  the  probability  distribution 
for  the  rband  source  size  is  log(ri  r /cm)  =  15. 3±0. 3,  where  the 
error  bar  represents  the  bounds  of  the  68%  confidence  interval. 
In  Figure  6  we  also  show  for  comparison  the  lcr  range  of  values 
for  the  accretion  disk  size  obtained  by  Mediavilla  et  al.  (2011, 
hereafter  Mil)  from  their  chromatic  microlensing  analysis  of 
SBS  0909  under  the  assumption  of  a  logarithmic  size  prior  and 
a  disk  temperature  profile  power-law  index  p  =  I  /a  =  4/3. 
In  order  to  make  the  most  accurate  comparison,  we  converted 
Ml  1  ’s  half-light  radius  to  a  thin-disk  scale  size,  scaled  the  result 
from  its  rest-frame  wavelength  of  1460  A  to  /.rcst(r)  =  2620  A 
assuming  R,  oc  /,4'/:1  for  thin  disks,  scaled  once  more  to  a  mean 
microlens  mass  (M)  =  0.3  M0,  and  corrected  for  inclination 
assuming  the  same  i  =  60°  we  applied  to  our  disk  sizes. 


As  can  be  seen  in  Figure  6,  our  microlensing  source  size  for 
SBS  0909  is  smaller,  but  marginally  consistent  with  Mil’s. 
We  suspect  that  the  difference  may  arise  from  the  effect  of  the 
magnification  pattern  pixel  sizes  on  the  size  distribution  (noting 
that  our  r  band  disk  size  is  similar  to  the  size  of  the  pixels 
in  Mil’s  magnification  patterns),  evidence  of  microlensing  of 
the  quasar’s  broad  emission  lines,  and  the  different  modeling 
approaches  used  in  the  different  studies.  It  is  also  conceivable 
that  Ml  l’s  result  has  been  affected  by  their  use  of  single-epoch 
spectra  and/or  their  combination  of  optical  and  near-IR  spectra 
obtained  at  epochs  separated  by  several  years:  the  uncorrected 
time  delay  and  intrinsic  variability  may  alter  the  continuum  and 
emission  line  flux  ratios  from  their  true  values. 

Because  it  is  possible  that  the  /-band  flux  which  we  ob¬ 
serve  from  SBS  0909  and  model  in  our  simulations  could  be 
contaminated  by  UV  or  optical  photons  scattered  by  the  broad 
line  region,  or  higher  energy  continuum  photons  reprocessed 
by  the  broad  line  region  and  reemitted  as  emission  lines,  our 
r-band  accretion  disk  size  may  be  an  overestimate  (see  Morgan 
et  al.  2010;  Guerras  et  al.  2013).  In  fact,  the  prominent  Mgn 
emission  line  (rest-frame  2798  A)  in  the  spectrum  of  SBS  0909 
falls  within  the  passband  of  our  r-band  filter  (see  Lubin  et  al. 
2000;  Mil),  so  contamination  from  line  emission  is  of  partic¬ 
ular  concern.  To  investigate  the  possibility,  we  have  repeated 
our  microlensing  simulations  under  the  assumption  that  a  frac¬ 
tion  of  the  observed  r- band  flux  should  actually  be  attributed  to 
unmicrolensed  emission  from  scattered  light  or  the  broad  line 
region.  We  find,  however,  that  unmicrolensed  contamination  is 
not  a  significant  factor  in  our  accretion  disk  size  determination: 
even  when  we  assume  that  as  much  as  30%  of  the  observed 
r-band  flux  is  contributed  by  contamination  from  emission  on 
large  physical  scales,  the  median  of  the  /  -band  size  probability 
distribution  is  essentially  unchanged  at  log(rSir /cm)  =  15.3/.0  4. 

We  show  the  posterior  probability  density  for  the  observed- 
frame  g-band  accretion  disk  size  for  SBS  0909  resulting  from 
our  two-band  microlensing  analysis  in  the  lower  panel  of 
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Figure  6.  Joint  probability  distributions  for  physical  scale  size  of  the  accretion  disk  in  observed-frame  r  (top  panel)  and  g  (bottom  panel)  bands  in  SBS  0909.  Both 
distributions  have  been  corrected  for  inclination  assuming  i  =  60°.  The  solid  and  dashed  vertical  lines  indicate  the  Schwarzschild  radius  and  the  radius  of  the  last 
stable  orbit  in  the  Schwarzschild  metric,  respectively,  for  black  holes  of  mass  10s  51  M@  and  109  29  Mq.  Our  disk  sizes  are  more  consistent  with  a  central  black  hole 
mass  of  108'51  Mq  for  SBS  0909.  For  comparison,  we  show  the  I  rr  range  for  the  disk  size  of  SBS  0909  obtained  by  Mediavilla  et  al.  (201 1),  scaled  to  our  rest-frame 
wavelengths  and  for  mean  microlens  mass  ( M )  =  0.3  Mq.  Our  disk  scale  radii  are  marginally  consistent  with  Ml  l’s  result. 


Figure  6.  We  have  included  in  this  panel  the  la  range  of  the  ac¬ 
cretion  disk  size  result  from  Ml  1  as  well,  corrected  as  described 
in  the  paragraph  above,  except  the  result  has  now  been  scaled 
to  the  rest- frame  wavelength  of  the  g  band  (2020  A)  instead  of 
r  band.  As  might  be  expected  due  to  the  shorter  time  baseline 
of  the  g-band  monitoring  data  and  somewhat  poorer  data  qual¬ 
ity,  the  constraints  we  obtain  for  the  g  band  size  are  not  nearly 
as  tight  as  those  for  the  r  band.  In  light  of  the  wide  peak  of  the 
g-band  probability  distribution  and  the  hints  of  secondary  peaks, 
we  regard  the  g-band  result  as  a  preliminary  estimate.  Despite 
the  large  uncertainty,  it  is  encouraging  that  the  median  of  the 
g-band  physical  size  distribution,  log(ri  g/cm)  =  14.8  ±  0.9,  is 
indeed  smaller  than  the  median  of  our  /  -band  size  distribution, 
consistent  with  the  shorter  wavelength  of  g-band  and  the  values 
for  the  disk  temperature  slope  in  the  literature.  Like  our  r  band 
disk  size,  our  g  band  result  is  notably  smaller  than  the  scaled 
result  from  Mil,  although  the  significance  of  the  discrepancy  is 
low  due  to  the  large  uncertainties.  We  also  calculate  the  proba¬ 
bility  density  for  the  ratio  of  the  r  and  g  band  disk  sizes,  which 
we  show  in  Figure  7.  The  distribution  is  very  broad,  with  a  me¬ 
dian  value  and  la  confidence  level  of  log  rsr/rSig  =  0.5^'q- 
Using  a  standard  thin-disk  temperature  profile  ( T  oc  A’  ,/4) 
to  predict  the  observed-frame  r/g  size  ratio  would  result  in 
log  rsr/rsg  —  0.15,  which  is  smaller  but  statistically  consistent 


with  our  observed  value.  Unfortunately,  our  best  value  for  the 
r/g  size  ratio  is  too  uncertain  to  provide  any  conclusive 
indication  of  the  temperature  profile  of  the  accretion  disk 
in  SBS  0909.  We  expect  that  a  future  analysis  utilizing 
a  longer  g-band  time  baseline  for  photometric  monitoring 
will  significantly  improve  the  precision  of  the  g-band  size 
measurement. 

An  interesting  point  made  by  Figure  6  is  that  the  accretion 
disk  sizes  we  obtain  from  our  two-band  microlensing  analysis 
are  not  consistent  with  the  central  supermassive  black  hole 
mass  derived  from  SBS  0909’s  H/i  emission  lines  by  Assef 
et  al.  (2011),  109-29  M0.  The  H/l  line  measurement  appears 
to  be  the  most  reliable  for  SBS  0909  in  Assef  et  al.  (2011), 
because  the  line  profiles  for  C  iv  and  Ha  emission  lines  were 
difficult  to  model.  By  extension,  we  might  expect  the  black 
hole  mass  calculation  from  H/i  to  be  the  most  reliable  as  well. 
However,  the  innermost  stable  orbit  of  a  maximally  rotating 
Kerr  black  hole  and  the  Schwarzschild  radius  predicted  for  a 
black  hole  of  mass  10929  M0  both  fall  within  the  la  bounds 
of  our  r-  and  g-band  accretion  disk  sizes.  Another  possible 
clue  that  the  H/l  black  hole  mass  may  be  problematic  is  that 
the  theoretical  scale  radius  of  a  thin  accretion  disk  at  2620  A 
surrounding  a  black  hole  of  mass  I  ()9  29  M0  (the  “theory  size,” 
log [ 7U620 /cm]  =  15.73)  is  larger  than  the  r-band  microlensing 
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Figure  7.  Joint  probability  density  for  the  ratio  of  the  accretion  disk  sizes 
in  observed-frame  r-band  and  g-band  (rs  r/rStg)  for  SBS  0909.  The  vertical 
line  highlights  the  location  of  rs  r/rs  g  =  1.  The  distribution  is  very  wide, 
reflecting  the  poor  constraints  we  are  able  to  place  on  the  observed-frame  g-band 
accretion  disk  size.  The  median  and  l<r  values  for  the  size  ratio  distribution  are 
log  rsr/rsg  =  0.5*°]q,  which  are  larger  but  not  statistically  inconsistent  with 
the  r/g- band  size  ratio  expected  for  a  thin  accretion  disk. 

size  from  our  simulations;  yet,  in  the  quasar  microlensing 
literature,  the  theory  size  is  consistently  smaller  than  the  results 
of  microlensing  simulations  (see,  e.g.,  Morgan  et  al.  2010; 
Blackburne  et  al.  2011).  The  accretion  disk  sizes  predicted  by 
our  microlensing  simulations  are  more  consistent  with  the  C  iv 
black  hole  mass  estimate,  10851  MQ,  highlighting  the  difficulties 
and  uncertainties  associated  with  estimating  black  hole  masses 
from  quasar  emission  lines. 

Our  Monte  Carlo  microlensing  analysis  of  SBS  0909  g-  and 
r-band  monitoring  data  has  enabled  us  to  estimate  the  r-to-g- 
band  accretion  disk  size  ratio  and  improve  the  precision  of  the 
system's  time  delay,  despite  uncertainties  in  the  macroscopic 
lens  galaxy  model.  We  suggest  that  deep,  high-resolution  imag¬ 
ing  of  the  SBS  0909  system  will  be  necessary  to  completely 
resolve  the  lingering  questions  about  the  lens  galaxy  model.  For 
the  moment,  our  delay  measurement  leads  us  to  favor  the  lens 
galaxy  astrometry  of  Lehar  et  al.  (2000)  over  that  of  Sluse  et  al. 
(2012).  Our  result  for  the  r/g- band  size  ratio  is  rather  coarse  and 
requires  confirmation;  to  that  end,  we  have  begun  monitoring 
SBS  0909  in  g-band  again.  Additionally,  we  recently  expanded 
the  USNA/USNO  lensed  quasar  monitoring  campaign  to 
near-IR  wavelengths.  We  look  forward  to  the  improvements 
in  our  ability  to  constrain  the  temperature  profiles  of  quasar  ac¬ 
cretion  disks  that  will  be  enabled  by  size  measurements  across 
a  significantly  larger  wavelength  baseline. 
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